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I. INTRODUCTION 



In quantum mechanics the exactly solvable models (ESM) are essential: since the com- 
plete physical information is encoded in few analytic expressions, they are ideal to test 
the convergence of numerical methods. Moreover, the ESM constitute the starting point 
for applying perturbation techniques. There are some other models, intermediate between 
the exactly solvable ones and those which can be solved just numerically, which are known 
nowadays as quasi-exactly solvable (QES). For them there exist analytical expressions en- 
coding just partial physical information, e.g., for only a part of the Hamiltonian spectrum 
Along the years several systems have been identified as belonging to the QES class, 
in particular, there was a dominant conviction that the associated Lame potentials were 
QES [2]. However, very recently it was shown that the associated Lame potentials for (m, t) 
integers are exactly solvable, in the sense that the stationary Schrodinger equation admits 
analytic solutions for any value of the energy parameter j^, 4]. The initial motivation to look 
for this result was the need to implement supersymmetric quantum mechanics for generating 
new exactly solvable models (periodic and asymptotically periodic) 0,0, 0, y,Q, Q. In 
order to implement non-singular transformations of general type, the explicit expressions for 
the Schrodinger solutions in the gaps as well as the band edge eigenfunctions were required 



In the next section we will discuss two techniques to find the general solution for the asso- 
ciated Lame equation for an arbitrary value of the energy parameter: first an ansatz based 
procedure and then a systematic technique related to the well known Frobenius method. 
In section 3 we will apply the supersymmetry transformations for generating new periodic 
and asymptotically periodic potentials which are almost isospectral to the initial associated 
Lame potential. In section 4 our general results will be illustrated through the particular 
case characterized by {m,i) = (3,2). Our conclusions will be finally given at section 5. 



II. GENERAL SOLUTIONS OF THE ASSOCIATED LAME EQUATION 



We want to find the general solution of the Schrodinger equation for the associated Lame 
potentials with an arbitrary value of the energy parameter E: 



H^{x) = [-dl + V{x)] ipix) = E'4j{x), 



V{x) = m{m + l)rsn^x + + l)k 



2 cn^x 
dn^x' 



(1) 
(2) 



where m, £ G N, m > i, and snx = sn(x, k), cnx = cn(x, k), dnx = dn(x, k) are the Jacobi 



elliptic functions of real periods AK, AK and 2K respectively 
Lame equation, which can be also expressed in Weierstrass form: 
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p{z) - ei 
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through the changes 



X — iK' _ ~ / \ r 

z = , Cj = ei — Cj, i = 2, 3, E = e^rn^m + 1) + [E 



l)]e3, 



p[z) = p{z, u, oj') being the Weierstrass elliptic function of half-periods u) = K/ y/e^, ui' = 

In order to solve ([3]), let us denote by 'ip'^{z) two linearly independent solutions. Thus, 
their product '^{z) = ip^ip~ will satisfy the following third-order equation: 
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p\z)m = 0. 



(4) 
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In the first place let us consider the following ansatz for 

A2 



(5) 



By plugging this in equation (jlj), the parameters Ai, A2, i,m can be fitted 

(a) m = 1, i=l, Ai = E + ei, ^2 = 6263, 
(6) m = l, £ = 0, Ai = E + ei, ^2 = 0. 

Once ^1/(2;) has been gotten, it is straightforward to find ip^{z) through the formula: 

In the non-trivial case (a) with (m, i) = (1, 1), in which the associated Lame equation is not 
directly the Lame one, the corresponding solution reads: 
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x-iK' 
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± a. 



„ , X-iK' I , , \ „ x-iK' 

a I — 7w=- + uJi) a ' 



exp 



where cui = uj, C'i.^) = ~p{^)^ l^^'^i^)]' = C(^)- 
On the other hand, for a modified ansatz of kind 

^{z) = [p{z) - ei]2 + B,[p{z) - ei] + + 



p{z) - ex 



it turns out that 
(c) m = 2, 
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- ei LBi + 6263, ^3 = 0. 



Once again, in the case (c) with {m,i) = (2, 1), in which the associated Lame equation is 
not the Lame one, the corresponding solutions are: 
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x-iK' 



±br 



„ , x—iK' I , , \ „2 / x~iK' 



exp 



Although the previous approach allows to find the solutions of equations ([3lll]) for some 
integer values of the parameters (m, £), however it is not completely systematic. In order to 



fill the gap, let us use the Frobenius method for solving (jl]) {4]. With this aim, let us make 
the following changes: 

ei - p{z) 



y 



62 



Therefore, the following equation for $ is obtained: 



(7) 



(8) 



where Pi{y) are the following i-th degree polynomials in y: 

P,{y) = 2y\e2y' -3e,y + 63), PM = 3y[e2{3-2i)y^ -Qe,{l-i)y +63(1-2^)], 

P^(y) = 2{e2(3£2 - - 6£ - m + 3)y^ - [E + ei{9f - - 9£ - m + 3)]y + £{2£ - 1)63}, 

Piiy) = 62(2^ -l){m + i){m-i + l)y + 2i[E + ei{3f - - m)]. 



Let us propose now 



from which the next equations arise 



r+p 



(9) 



r=0 



ao/o(p) = 0, ai/o(p + 1) + ao/i(p) = 0, 
ar+2fo{p + r + 2) + ar+ifi{p + r + 1) + ar/2(p + r) = 0, 

/o(p) = e3p(p - 1 - 2£)(2p -2i- 1), /i(p) = 2{p - £){ei[m(m + 1) - 3(p - 
/2(p) = e2(p -m-e){p + m-£ + l){2p+l- 2i). 



From the three roots p = 0, p = 2£ + l, p = £+ l/2of the indicial equation /o(p) = 0, the 
series Qj can be made finite just for p = 0, since for r = 2£ — 1 it turns out that 



a2,/i(2£) + a2^_i/2(2£-l) = 



where the determinant family 



-ao, r= 1,2,...,2£, (10) 



n Ms) 

s=l 



/i(r-l) /2(r-2) 
/o(r-l) /i(r-2) /2(r-3) 
/o(r-2) /i(r-3) /2(r - 4) 







/i(0) 



(11) 



satisfies the following recurrence relationship: 

F, = /i(r-l)F,,_i-/2(r-2)/o(r-l)F,_2, r G N, Fo = l, F_,=0. (12) 
In addition, we have to fix a2i+i so that the series ends up after (A^ + 1) terms, i.e., 

a^/i(iV) + a^_i/2(iV - 1) = 0, ajv+i = 0, ^(iV) = 0. 
The previous equations imply that 

N = m + i, am+efiim + £) + am+i-if2{m + i-l) =0, am+e+i = 0. 
Two different cases can be identified: 

(i) For m = £ it turns out that formula (fTOj) provides all non-null a^, r = 0, 1, . . . , 2£. 

(ii) For m = i + u, u = 1,2, . . ., equation (ITOl) supplies just Qr, r = 1, . . . , 2£, while the 
remaining non-null coefficients are given by: 



a2£+r = a2£, r = 1, 2, ... I/, 

where is the minor of -Pj^+^+i-r in Laplace expansion of the determinant F2e+u+i- 
By coming back now to equation ([7]), up to a constant factor in both cases we have 



(13) 



m+e 

n MZ) - pibr)] 
r=l 



m+e 

p{bi), p{b2), . . . , p{bm+i) being the zeros of (^Ai^i —t)/(^2Y- By using once again equation 

r=0 

IQ, the solutions of the associated Lame equation are finally: 
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exp 
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(14) 



Note that the well known solutions of the Lame equation are recovered for £ = 0: 
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III. SUPERSYMMETRY TRANSFORMATIONS 

The supersymmetric quantum mechanics, as an approach to generate new exactly solvable 
Hamiltonians H from an initial solvable one H, is based on the intertwining relation 

HB = BH. (15) 

This means that the eigenf unctions of H are constructed through the non-null action of B 
onto the eigenfunctions of H . If 5 is a fc-th order differential operator it turns out that 
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V{x) = V{x) - 2[lnPF(ui, . . . (16) 

where Ui are solutions, which can be nonphysical, of the stationary Schrodinger equation 
associated to k different factorization energies ei,i = 1, . . . , k, i.e. 

Hui = eiUi. (17) 

For periodic potentials, it has been realized that when Bloch-type seed solutions associated 
to factorization energies which belong to the energy gaps are used, then the SUSY partner 
potentials V{x) are again periodic, isospectral to the initial one Baa 

. On the other hand, 

when appropriate linear combinations of the solutions are employed, the SUSY partner 
potentials V{x) of V{x) become asymptotically periodic, with periodicity defects appearing 
due to the creation for H of bound states embedded into the gaps j^, 9]. This suggests 
a natural ordering for the SUSY transformations of periodic potentials which will be next 
followed (we restrict ourselves to first and second-order techniques). 



A. Periodic first-order SUSY partner potentials 

Let us take u{x) = ip'^i^x) (see dHl)) with e < Eq, where Eq is the lowest band edge 
eigenvalue for the corresponding associated Lame potential. The periodic first-order SUSY 
partner potentials of V{x) thus read 
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B. Asymptotically periodic first-order SUSY partner potentials 



Let us choose now u{x) as a general linear combination of the solutions ip^ (e < £"0) 



u{x) = Aip^{x) + Bip (x) = Aip^ {x)(t)^{x) = Bip (a;)0_(a;), 



a 



■x-iK' 



K ^ A. 



(19) 



r=l " ^ V^S 

The first-order SUSY partners for the associated Lame potentials become now asymptotically 
periodic, with explicit expressions given by 



y"P(a;) = V±(x) - 2[ln0±(x)]", 



(21) 



where V±{x) is given by f|T8|) . The spectrum of the Hamiltonian iJ"^ contains the allowed 
energy bands of H but in addition it has an isolated bound state at £^ = e. 

C. Periodic second-order SUSY partner potentials 

Let us take now two factorization energies ei^2 inside the same energy gap, and the 
corresponding Schrodinger seed solutions in the way 

X 



n + K 

r=l 



U\{x) = i)'! ix) = „, -r^i , , -wi 

ivy ) 



exp 



ii2(a;) = ■?/'^(x) 



n 

r=l 



r=l 



a^(^ + ^0^-(^) 



exp^[£C(^i)-EC(^- 



r=l 



It turns out that the Wronskian of M1.2 is nodeless, which is conveniently expressed as: 



W (mi, M2) = V'J*" (a;) '^t {x) g{x) g{x) 



In 



The second-order SUSY partner potentials of V{x) become periodic: 



V{x) = m{m - 3)khn^x + i{i - S)k'^^ 



m+l 

+2^ E [sn2(x + v^6,) + sn2(x + y^^K)] - 2 (In g)" . (22) 

r=l 
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D. Asymptotically periodic second-order SUSY partner potentials 



Let us choose once again ei^2 in the same energy gap but now -Ui 2 are general hnear 
combinations of 



2 ; 



where 



^1 



^/ x—iK' 
r=l ^ ves 

m+i ^( x-iK 



+ 6J 



■x-iK' 



r=l y/e3 

An appropriate choice of A^2 leads to a nodeless Wronskian, which is expressed as: 



W{Ui,U2) = MlM2fl'"^ = 1pt'^2(Pl(p2 9 



np 



Ul 



In 



+ 



In 



9 + 



In 



The second-order SUSY partner potentials are again asymptotically periodic: 



= V{x)-2 

where V{x) is given by fl22l) . 



In 



'1 V2 



(23) 



(24) 
(25) 



(26) 



IV. EXAMPLE 



In some previous papers it has been studied the associated Lame potential 



partners for (m,£) = (1, 1), (m,£) = (2, 1), and (m, = (3, 1) Q, 4, [isl, [jj]. Here we will 
illustrate our general procedure for the associated Lame potentials with (m, = (3,2), i.e.. 



s and its SUSY 



V{x) = 12khn^x + 6A;2^. 

dn X 



(27) 



Note that there are explicit expressions for the band-edge eigenfunctions and eigenvalues of 



H 



14| : in particular for the 'ground' state it turns out that: 



^o(^) = dn^^, Eq = 9P. 



(28) 
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Since m + £ = 5, we have to evaluate five constants 01,02,03,04,05, (without losing 
generality we have taken oq = 1). Let us write down first the basic elements /j. 



Up) = espip - 5)(2p - 5), = 2(p - 2)[3eip(4 - p) - E], 

f2{p) = e2{p - 5)(p + 2)(2p - 3), E = e,{E - 6) + Ue^. 

Note that the four constants Oi, 02, 03, 04 are determined from ffTOl) 



Oi 



02 



/o(l)' /o(l)/o(2 
while 05 is computed from ( |T3T 



03 



05 



/o(l)/o(2)/o(3) 
^^0/2(4) 



, 04 



/o(l)/o(2)/o(3)/o(4)^ 



-04. 



(29) 
(30) 

(31) 
(32) 



From film 121) one may calculate quite straightforwardly the Fj's 

Fi = 4E, F2 = 8(^2-961^-456263), F3 = -2886263^, F4 = 5(726263)'. (33) 



To obtain the D,'s we need 



/i(5) /2(4) 
/o(5) /i(4) /2(3) 















/o(4) /i(3) /2(2) 

/o(3) /i(2) /2(1) 

/o(2) /i(l) /2(0) 

/o(l) /i(0) 



(34) 



Then 



Dq = [minor of Fq in Fq] = 1, 
Di = [minor of F5 in Fq] = /i(5) 

Finally, using the recurrence relation ([T^ we obtain from 



(35) 
(36) 



E E^ - 9eiE - 456263 

flO = 1, fll = -77—5 ^2 = 

363 



9e1 



03 



262^ 
' 361 ' 



5e2 



25e^ 



^2 ' 



05 



ei(E + 156i) 



(37) 



We employ these coefficients to find the roots c^, r = 1, . . . , 5 of the fifth-order equation 

5 



r=0 



61 — t 
62 



0. 



(38) 
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These roots are used then to invert the transcendental equation p{br) = Cr to determine the 
6r's (with the restriction '$'\z=br > 0), which are thus inserted in the exphcit expressions 
for ip'^lx). Finally, the resulting Bloch solutions can be used, either directly or in the 
corresponding Wronskian, to derive the periodic SUSY partner potentials V±{x) of ( |T8i) or 
( l22ll . On the other hand, different linear combinations of kind (fT9l) or (l23ll can be used 
to derive the potentials V"'^{x) of (!2T!) or (!26|) . which have periodicity defects. The final 
results of these procedures are illustrated in figures below, in which we show in gray the 
original associated Lame potential for m = 3, £ = 2, = 0.9. In Fi gure la we plot in 
black one of its periodic first-order SUSY partners generated through a Bloch solution with 
e = 8<-Eo = 8.1, while in Figure lb it is illustrated one of its asymptotically periodic 
partners for the same e. On the other hand, in Figure 2 we have drawn similar graphs 
(black curves) for the corresponding second-order SUSY partner potentials, periodic and 
asymptotically periodic. For the periodic case (Figure 2a) we have used two Bloch solutions 
Ui{x) = ipilx), U2{x) = ipti^) associated to ei = 10 and €2 = 10.1, which fall in the first 
finite energy gap (8.1031, 11.7154). For the asymptotically periodic case (Figure 2b) we have 
used the same pair of factorization energies, with linear combinations Ui{x) = '?/'J*"(x) + '?/'f (x) 
and U2{x) = 1^2 i^) "~ 1-5 4^2 {x). Note that in both asymptotically periodic cases, of first 
and second order, the periodicity defects are clearly detected. 
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Figure 1: First-order SUSY partners (black curves) for the associated Lame potential (gray curves) 
with m = 3, i = 2, k"^ = 0.9. (a) Periodic case generated through the Bloch solution ■0"^(x) with 
e = 8. (b) Asymptotically periodic case in which the linear combination ■i/;"^(x) -|- ip~ (x) for the 
same e is used. 
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V. CONCLUSIONS 



We have shown that the associated Lame potentials for integers values of the parameter 
pair (m, i) are exactly solvable. When using as seeds Bloch-type solutions inside the gaps 
new exactly solvable periodic potentials are generated. On the other hand, for seeds chosen 
as general linear combinations of Bloch-type solutions, the SUSY partner potentials become 
asymptotically periodic, the corresponding spectra having bound states embedded into the 
gaps. These potentials are interesting since the new levels could work as intermediate tran- 
sition energies for the electrons to jump between the energy bands. 

(a) (b) 





Figure 2: Second-order SUSY partners (black curves) for the associated Lame potential (gray 
curves i with m = 3, e = 2, = 0.9. (a) Periodic case generated through two Bloch solutions 
ui{x) = ipf{x), U2{x) = tp2ix) for ei = 10, 62 = 10.1. (b) Asymptotically periodic case in which 
the linear combinations ui{x) = 'ipf{x) + ip^ (x), U2{x) = ip2i^) ~ ^-^ "4^2 i^) same ei^2 are 

used. 



We would like to end up this paper by making a historical precision concerning SUSY 
techniques applied to periodic potentials (see also the discussion in [15|)- It is a fact that the 
recent interest on the subject [3, Z S, I^, Q, I^, 10, 1^, 17, 18] was catalyzed by Dunne and 
Feinberg discovery of the self-isospectrality for the Lame potentials with m = 1, induced by 
the first-order SUSY transformation which employs as seed the ground-state eigenfunction. 
Selt-.ospect.al,ty .neans, in particular, that the SUSY pattnet potent.al beco„.es just a 
displaced version of the initial one, by half the period in [5]. Soon it was realized that the 
self-isospectrality in which the new potential becomes the initial one displaced by any real 
number arises as well for the Lame potential with m = 1, the seeds employed being Bloch- 

n Q 

type solutions associated to factorization energies in the gaps [8|, 19|. More general SUSY 
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transformations, either of higher order or involving general solutions of the Schrodinger 



equation for a given factorization energy, have been also introduced 



However, it 



is worth to note that there are several interesting works, previous to Dunne and Feinberg 
paper, in which the SUSY techniques were applied to periodic potentials (see e.g. In 
particular, it is remarkable the work of Braden and Macfarlane in which the self-isospectrality 
for the Lame potential with m = 1 was discovered for the first time 19|. This is a typical 
story of a discovery followed by a later rediscovery, which often arises in science. Our opinion 
is that both works are valuable, complementary to each other, and hence they are worth to 
be studied in detail. 
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